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Abstract 

The search for supersymmetry at Run 1 of LHC has resulted in gluino mass limits nig >1.3 TeV 
for the case where mq 3> nig and in models with gaugino mass unification. The increased energy 
and ultimately luminosity of LHC13 will explore the range mg ~ 1.3 — 2 TeV. We examine how 
the discovery of SUSY via gluino pair production would unfold via a comparative analysis of three 
LSP archetype scenarios: 1. mSUGRA/CMSSM model with a bino-like LSP, 2. charged SUSY 
breaking (CSB) with a wino-like LSP and 3. SUSY with radiatively-driven naturalness (RNS) 
and a higgsino-like LSP. In all three cases we expect heavy-to-very-heavy squarks as suggested by 
a decoupling solution to the SUSY flavor and CP problems and by the gravitino problem. For 
all cases, initial SUSY discovery would likely occur in the multi-6-jet + fdj, channel. The CSB 
scenario would be revealed by the presence of highly-ionizing, terminating tracks from quasi-stable 
charginos. As further data accrue, the RNS scenario with 100-200 GeV higgsino-like LSPs would 
be revealed by the build-up of a mass edge/bump in the OS/SF dilepton invariant mass which 
is bounded by the neutralino mass difference. The mSUGRA/CMSSM archetype would contain 
neither of these features but would be revealed by a buildup of the usual multi-lepton cascade 
decay signatures. 
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I. INTRODUCTION 


The LHC8 (LHC with a/s = 7 — 8 TeV) era has come to a close and the LHC13 era is 
underway! What have we learned from LHC8? The Standard Model (SM) has been spec¬ 
tacularly confirmed in a vast assortment of measurements [Tj. And most importantly, a very 
SM-like Higgs boson has been revealed with mass rrih = 125.09 ± 0.24 GeV (ATLAS/CMS 
combined) (2j [3j . The next major target for LHC is to root out evidence for supersymmetry 
(SUSY). Indeed, it has been declared that if LHC13 does not uncover evidence for SUSY 
early within Run 2, then physics will have entered a state of crisis [3]! 

What we have learned from LHC8 is that- in generic models such as mSUGRA/CMSSM- 
no evidence for SUSY translates into mass bounds of 

rrig > 1.3 TeV for rrig and (1) 

rrig > 1.8 TeV for rriq ~ rrig. (2) 

In addition, the rather large value of rrih — 125 GeV seems to require large radiative cor¬ 
rections to m 2 in the MSSMj5j. The Higgs mass can be accommodated with TeV-scale top 
squarks for large trilinear soft breaking parameter Ao[6], or by 10-100 TeV top squarks in 
the minimal mixing case[7j. Naively, these rather high sparticle mass limits seem to conflict 
with notions of weak scale naturalness which favor sparticles at or around m wea k — 100 GeV, 
the value of m-wz,h- This has led to some puzzlement as to the emerging Little Hierarchy: 
why is m(sparticle) ^ ■ m wea k ? It has also led to more detailed examination of what is meant 
by electroweak naturalness. 

The point of contact between SUSY Lagrangian mass parameters (soft terms and super¬ 
potential n term) and hard data occurs in the scalar (Higgs) potential: in the MSSM, it is 
given by 

Vffiggs Uree + AV, (3) 

where the tree level portion is given by 

Vtree = + P) Kt + « + S) l&SP 

-BuKK + h.c.) + t(g 2 + g J2 )(K I 2 - |ft»| 2 ) 2 (4) 

and the radiative corrections (in the effective potential approximation) by 

Al/ = E 

where the sum over i runs over all fields that couple to Higgs fields, is the Higgs field de¬ 
pendent mass squared matrix (defined as the second derivative of the tree level Lagrangian) 
of each of these fields, and the trace is over the internal as well as any spin indices. Mini¬ 
mization of the scalar potential in the h° u and h Q d directions allows one to compute the gauge 
boson masses in terms of the Higgs field vacuum expectation values v u and Vd, and leads to 
the well-known condition that 
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where the E“ and E([ terms arise from derivatives of AV evaluated at the potential minimum 
and tan B = —. This minimization condition relates the Z-boson mass scale to the soft SUSY 

r Vd 

breaking terms and the superpotential higgsino mass p. In most computations of the SUSY 
mass spectrum, the weak scale soft terms are determined by renormalization group running 
from a constrained set of parameters set at some high scale A. In gravity-mediation[8], A 
is usually taken to be maur — 2 x 10 16 GeV[9]. Then the weak scale value of p is dialed 
(fine-tuned) so that the measured value of mz is obtained. An evaluation of the extent of 
this fine-tuning is provided by the electroweak measure A ew which evaluates the largest of 
the 43 terms on the right-hand-side of Eq. [6j If one term on the RHS is m 2 z , then some 
other unrelated term will have to be large and of opposite-sign to guarantee that rri z = 91.2 
GeV. To avoid such large weak scale tuning, evidently all terms on the right-hand-side of 
Eq. [6] should be comparable to or less than m|. This implies the following]!]], TT] , 

• The superpotential higgsino mass p ~ 100 — 200 GeV, the closer to mz the better. 
The lower limit p > 100 GeV comes from null searches for chargino pair production 
at LEP2. 

• The soft term m? Hu is radiatively driven to small values ~ at the weak scale. 

• The radiative corrections E“ are not too large. The largest of these usually comes 
from the top-squarks. Each of the terms E^U^) are minimized by TeV-scale highly 
mixed top squarks, a condition which also lifts mh up to ~ 125 GeV[[10], fTl ] . 

Some alternative fine-tuning measures also have been advocated in the literature. 

1. The usual application of the Higgs mass large-log measure A hs = dm 2 H / {m\/ 2) where 
5m 2 Hu ~ Aj ( m Q 3 + m u 3 + A 2 ) In (A /msusy) has been challenged [T2l ITS] in that it ig¬ 
nores the dependent term m 2 H which occurs in the RGE. But the larger m 2 H { A) 
becomes, the greater is the cancelling correction to 8mjj [TT]. By appropriately com¬ 
bining dependent terms, A hs reduces to the same general consequences as A ew- 

2. Alternatively, the Ellis et a/./Barbieri-Giudice measure[151 [16] is defined as A bg = 
max,; \d ln(m 2 z ) / dpi\ where the p, constitute fundamental high scale parameters of the 
theory. To evaluate A bg, ttt.| must be evaluated in terms of fundamental high scale 
parameters usually taken to be the GUT scale soft breaking terms. The usual ap¬ 
plication of this measure has been challenged[TT1 [T5] in that in supergravity theories, 
the soft terms are not independent, but are evaluated as multiples of the fundamental 
gravitino mass 7713 / 2 - Evaluating A bg m terms of the independent parameters p and 
777 . 3 / 2 , low A bg also leads to the same general consequences as Aew- 

Using A ew, then indeed most constrained high scale SUSY models are found to be highly 
tuned in the EW sector [133 . An exception occurs for a pocket of parameter space of the 
two-extra parameter non-universal Higgs models [TT] where /i ~ 100 — 200 GeV and where 
m 2 Hu is driven to small negative values comparable to —m 2 z while allowing for highly mixed 
TeV-scale top squarks which provide ~ 125 GeV. This pocket of parameter space we call 
SUSY with radiatively-driven naturalness, or RNS for short. By requiring EW naturalness, 
then upper bounds can be computed for all sparticle masses m- In radiative natural SUSY 
with A ew < 10 (30) then it is found that [TT]: 

• Trig < 2 - 5 ( 5 ) TeV ’ 
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• rn h < 2 (3) TeV, 

• m wi,z 1 , a % 200 (300) GeV. 

The first of these values can be compared to the ultimate reach of LHC14 with 1000 fb _1 
where a 5cr discovery can be established for rrtg < 2 'itevpsmg. Thus, while EW naturalness 
certainly allows for gluinos and squarks to lie well beyond the ultimate reach of LHC13, it 
is also true that the most natural values of gluino and squark masses are those within the 
exploratory range of LHC13: the lighter the better. This motivates an examination of how 
a SUSY discovery via gluino pair production is likely to unfold at LHC13 when the gluino 
mass lies just beyond present bounds. 

In this paper we assume a gluino mass of rrig = 1400 GeV, i.e. just beyond present bounds. 
We then investigate how a SUSY discovery would unfold under three lightest SUSY particle 
(LSP) archetype scenarios: 

• the mSUGRA/CMSSM model with a bino-like LSP, 

• a charged SUSY breaking (CSB) scenario with a wino-likc LSP and 

• SUSY with radiatively-driven naturalness and a higgsino-like LSP. 


Our goal is to look for commonalities and differences between these three archetype scenarios 
that would allow a rapid determination of the nature of the LSP if a gluino pair production 
signal emerges at LHC13. 

Towards this end, in Sec. [HJ we present three archetype benchmark models (BM) labelled 
as mSUGRA, CSB and RNS. While each BM model contains a gluino with mass 1400 GeV, 
their implications for collider searches will be very different. In Sec. Ill we discuss how 
SUSY discovery would unfold in each BM model while in Sec. [TV] we discuss how each 


archetype could ultimately be distinguished as more integrated luminosity accrues. Briefly, 
in all cases the most likely initial discovery channel could occur in the rnulti-6-jet + $ T 
channel with ~ 3—8 fb _1 of integrated luminosity. For the CSB benchmark, the model would 
be distinguished by the presence of one or more cm-length highly ionizing tracks (HITs) from 
quasi-stable charginos which are produced within the gluino cascade decays. For the RNS 
scenario, as 100-1000 fb _1 of integrated luminosity accumulates, then a distinctive opposite- 
sign/same flavor (OS/SF) dilepton invariant mass edge should develop in multi-6-jet + $ T 
events which contain such a dilepton pair. The mass edge occurs at the kinematic limit 
m(£ + £~) < nig — ~ 10 — 30 GeV in RNS models. For the mSUGRA benchmark, 


neither of the above distinctive features should develop, but instead the usual multi-lepton 
plus multijet + $ T cascade decay topologies should build up as greater integrated luminosity 
accrues. Our summary and conclusions are given in Sec. |V} 


II. BENCHMARK MODELS 

In this section, we present three benchmark models representing each of three LSP 
archetype scenarios. Each scenario contains a light Higgs scalar mh — 125 GeV 1 and a 


1 We allow for a ~ ±2 GeV theory uncertainty on the Isajet RG-improved one loop effective potential 
calculation of rrih which includes leading two-loop terms [5]. 
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gluino of mass rn-g = 1.4 TeV, just beyond the bounds from LHC8. All spectra were gener¬ 
ated using the Isajet/Isasugra 7.84 program[20]. 


A. mSUGRA/CMSSM 

In the minimal supergravity model (mSUGRA or CMSSM)[8], it is assumed that super¬ 
gravity is broken in a hidden sector leading to a massive gravitino characterized by mass 
777-3/2, with 7773/2 ~ 1 TeV in accord with phenomenological requirements. I 11 the limit as 
M P —> 00 but keeping 7773/2 fixed, then one is lead to the global SUSY Lagrangian of the 
MSSM augmented by soft SUSY breaking terms each of order 7773 / 2 - A simplifying assump¬ 
tion (with minimal motivation) is that all soft scalar masses are unified to mo at the GUT 
scale. In addition, all gaugino masses are unified to 7771 / 2 , all trilinears are unified to Aq 
and there is a bilinear term B. Renormalization group running connects the GUT scale 
parameters to the weak scale ones. At the weak scale, the scalar potential is minimized and 
the superpotential p parameter is dialed (fine-tuned) so as to generate the measured value 
of m z = 91.2 GeV. 

Spectra from this popular model [?T1 BTj can be generated with many computer codes. 
In Table |T we show a mSUGRA benchmark model with mo = 5 TeV, 7771/2 = 517 GeV, 
Aq = —8 TeV and the ratio of Higgs vevs tan/l = 10. These parameters lead to a spectra 
with a gluino mass mg = 1.4 TeV, i.e. just beyond the reach of LHC8. The light Higgs mass 
?77h = 123.6 GeV, in accord with its measured value if one allows for the ±2 GeV uncertainty 
in our calculation of ?77ft. The Z\ is a bino-like LSP. The superpotential p parameter turns 
out to be p = 2861 GeV leading to A EW = 1968 so that this benchmark is highly fine-tuned 
in the EW sector. The calculated thermal neutralino abundance Vi~ p h 2 = 317, far beyond 

the measured value. Thus, some sort of 1. late entropy dilution, 2. decay of Z\ to an 
even lighter LSP such as an axino or 3. A-parity violating decays of Z\ would need to be 
invoked to bring the model into accord with the measured dark matter density. A schematic 
illustration of the lighter spectral states of the mSUGRA benchmark is shown in Fig. [l] 


B. Charged SUSY breaking 

In models labeled as minimal anomaly-mediation (mAMSB)|25j. it is assumed that SUSY 
is broken in a secluded sector so that the dominant contributions to soft terms come not from 
tree-level supergravity but from the superconformal anomaly. Such models leads character¬ 
istically to spectra including wino-like gauginos as the lightest SUSY particles[26j. Further, 
one obtains spectra with well-known tachyonic sleptons. In the original construct [25], it 
was suggested to augment soft scalar masses with a common m q term to cure the tachyon 
problem. 

The original mAMSB models seem disfavored in that they have problems generating 777ft ~ 
125 GeV due to a rather small weak scale A t soft term)?], 2?]. An alternative incarnation 
goes under the label of PeV SUSY[SH], split SUSY[22], pure gravity mediation[50] and spread 
SUSY[31]. In the simple yet elegant construction of Wells[2H], it is argued that the PeV scale 
(with 777 (scalars) ~ 7773/2 ~ 1 PeV=1000 TeV) is motivated by considerations of wino dark 
matter and neutrino mass while providing a decoupling solution [32] to the SUSY flavor, CP, 
proton decay and gravitino/moduli problems. This model invoked “charged SUSY breaking” 


5 




mSUGRA 

CSB 

RNS 

mo 

5,000 

50,570 

5,000 

Mi 

517.0 

927.3 

517.8 

m 2 

517.0 

140.5 

517.8 

M 3 

517.0 

-421.5 

517.8 

Ao 

-8,000 

140.5 

-8,000 

tan / 3 

10 

10 

10 

h 

2,861 

2,000 

150 

m A 

5,666 

2,000 

2,000 

m h 

123.6 

126.4 

124.1 

nig 

1,400 

1,399 

1,399 

mu L 

5,065 

50,205 

5,038 

m h 

1,929 

34,327 

1,332 

m w 2 

2.872.0 

2,064.8 

464.3 

m w t 

460.8 

143.6 

150.7 

m z 4 

2,866.3 

2,062.8 

473.6 

m z 3 

2,865.1 

2,062.2 

243.3 

m z 2 

459.8 

438.8 

159.5 

m z ! 

234.3 

143.4 

132.1 

Bino frac. 

0.9999 

0.0022 

0.2915 

Wino frac. 

0.0010 

0.9993 

0.1747 

Higgsino frac. 

0.0151 

0.0365 

0.9405 

n T ~ p h 2 

317 

0.0013 

0.01 

&-EW 

1968 

5228 

10.4 


TABLE I: Input parameters and masses (in GeV) for three benchmark points computed with 
ISA JET 7.84 |20| . Also displayed are the bino, wino and Higgsino fractions. 


(CSB) where the hidden sector superfield S is charged under some unspecified symmetry. 
In such a case, the scalars gain masses via SUGRA 




4Fs 

Ml ' 


(7) 


while gaugino masses, usually obtained via gravity-mediation as 

[ d 2 e^-WW => -^AA, (8) 

JMpMp w 

are now forbidden. Then the dominant contribution to gaugino masses comes from AMSB: 


Mi 


33 g\ 

5 167T 2 


m 3/2 ~ m 3 / 2 /120, 


(9) 
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FIG. 1: Spectra of benchmark models: mSUGRA (left), CSB (middle), RNS (right). We display 
the bottom part of the spectra up to 4 TeV where the electroweak -ino’s are shown in blue while the 
gluinos, stops and sbottoms are shown in red. The sparticles within the same column are ordered 
in increasing mass from left to right. In the CSB model, the stop and sbottom masses are ~ 35 
TeV: see Table [0 


M 2 = Y^m 3/2 ~ m 3/2 /360, (10) 

M 3 = _3 ^V m 3/2-"ks/2/40. (11) 

Saturating the measured dark matter abundance with thermally-produced (TP) winos re¬ 
quires rriy-r ~ M 2 ~ 2.5 TeV which in turn requires the gravitino and scalar masses to occur 
at the ~ 1000 TeV (1 PeV) level. A virtue of the CSB model is that the highly massive top 
squarks m$ ~ 50 — 100 TeV lead to rrih ~ 125 GeV even with a tiny A t trilinear soft term. 

The CSB benchmark point is listed in Table IT] where mo — m 3 / 2 = 50.57 TeV leading 

to squark and slepton masses ~ 50 TeV but with nrig = 1.4 TeV. The LSP is a wino-like 

neutralino Z\ with mass mg = 143.4 GeV. The superpotential /i parameter is taken to be 
2 TeV. The dominant contribution to the EW fine-tuning measure /S.ew comes from the 
top squark radiative corrections leading to A ew — 5228 so the model is highly fine-tuned 
in the EW sector. The thermally produced wino-like neutralino abundance is found from 
IsaReD 1551 to be = 0.0013 so WIMPs are thermally underproduced. They could be 
augmented via non-thermal WIMP production ( e.g. from gravitino, axino, saxion or moduli 
decays [34]) or the DM abundance could be augmented by other species such as axions[35]. 
The CSB benchmark is also shown schematically in Fig. [l] 


C. SUSY with radiatively-driven naturalness (RNS) 

In models with radiatively-driven naturalness, it is assumed that soft terms arise via 
gravity mediation and are characterized by the scale m 3 / 2 ~ 2 — 20 TeV. Such heavy soft 
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terms lead to rrih ~ 125 GeV for highly mixed TeV-scale top squarks. The // parameter 
arises differently. In the SUSY DFSZ axion model[361 [37] . the Higgs multiplets H u and 
Hd are assigned PQ charges so that the usual // term is forbidden although now the Higgs 
superfields may couple to additional gauge singlets from the PQ sector. The /j term is 
then re-generated via PQ symmetry breaking at a value of /j rv_/ fl/Mp so that the Little 
Hierarchy fi <C 777 . 3/2 is merely a reflection of the mis-match between PQ breaking scale 
and hidden sector mass scale f a -C m hidden . I 11 the MSY SUSY axion model[55], the PQ 
symmetry is broken radiatively as a consequence of SUSY breaking in a similar manner that 
EW symmetry is radiatively broken as a consequence of SUSY breaking. The radiative PQ 
breaking generates a small // ~ 100 — 200 GeV (as required by naturalness) from multi-TeV 
values of 7773/2 [39j. Once /i is known, then the weak scale value of m 2 Hu is determined by the 
scalar potential minimization condition and is also of order — m| as required by naturalness. 
The weak scale value of m 2 H is evolved to m GUT where it is found that mH u { m GUT ) % 777 0 , 
where m^ now labels just the matter scalar masses. 

The RNS benchmark model is shown in Table [I] with matter scalar mass rn () = 5 TeV and 
a trilinear soft term A 0 = —8 TeV. The ratio of Higgs vevs tan /3 = 10 and the pseudoscalar 
Higgs mass mA is taken as 2 TeV. The unified gaugino mass mi /2 = 517.8 GeV leading 
to 777 g = 1.4 TeV. The highly mixed top squarks with mass 777^ = 1.3(3.5) TeV lead to 

777ft = 124.1 GeV. Since p = 150 GeV, then the model has A ew = 10.4 or about 10% 
fine-tuning in the EW sector: the model is very natural. The LSP is a higgsino-like WIMP 
with mass m% i = 132.1 GeV. The TP relic density H~ p /7 2 = 0.01 but in this case the axion 

could comprise the bulk of DM|40]. 2 The RNS benchmark is schematically shown as the 
third frame of Fig. [l] 

III. HOW SUSY DISCOVERY UNFOLDS 

A. Gluino pair production 

In the benchmark scenarios we have selected, a heavy spectrum of matter scalars- squarks 
and sleptons- is assumed. This is in accord with at least a partial decoupling solution to 
the SUSY flavor, CP, gravitino and proton-decay problems. In addition, to accommodate 
Affleck-Dine[42] leptogenesis, then a non-flat Kahler metric is required[l3j from which one 
would expect generic flavor and CP violation. The decoupling solution allows the AD mech¬ 
anism to proceed in the face of potential flavor violations. 

In the case of decoupled matter scalars, then we expect gluino pair production and 
possibly electroweak -ino pair production to offer the main SUSY discovery reactions. In 
Fig. [ 2 J we show the NLO values of a(pp —> ggX ) reaction versus m^ for i/s = 8, 13 and 
14 TeV. The squark masses have been set to 5 TeV. We use Prospino to calculate the total 
cross sections[44j. 

For our benchmark points with mg = 1.4 TeV, we see that the LHC8 total production 
cross section a(gg) is about 0.6 fb. As A /s is increased to 13 TeV for LHC Run 2, then the 
total gluino pair production cross section jumps by a factor of ~ 30 to ~ 20 fb. Future LHC 
runs with fully trained magnets may attain y/s ~ 14 TeV for which cr(gg) would rise to ~ 35 


2 An alternative way to match the measured DM density is to reduce the bino mass Mi for the case of 
gaugino mass non-universality: see m ■ 



fb. While EW -ino pair production rates should be comparable to gluino pair production- 
due to their lower masses- we expect at this stage that gluino pair production is more easily 
seen due to its large energy release and no cost for leptonic branching fractions in the major 
signal channel of jets + $ T . 



m- (GeV) 

FIG. 2: Gluino pair production cross section at the LHC for y / s = 8,13,14 TeV calculated at NLO 
with Prospino[43j. Squarks are assumed to be heavy with mass rriq = 5 TeV. The shaded areas 
show the scale uncertainty. 


B. Gluino branching fractions and signatures 

Once produced, the gluinos can cascade decay[3Ej to a variety of final states which are 
listed in Table [TI} The decay modes including q in the final state are summed over q = u,d,s,c 
possibilities. It is evident from the Table that in all cases the decays to third generation 
quarks are enhanced over Erst /second generation quarks. Gluino three-body decays to third 
generation quarks were first calculated in Ref’s [461 |48] where their enhancement was noticed 
to arise from 1. couplings which include the large 6 and t Yukawa couplings, 2. generically 
smaller mediator masses < rriq and 3. large L-R mixing effects. For our benchmark 
models, we see that in mSUGRA, the g decays to states including 66 (both directly and 
via decay to top followed by t —> bW ) 81% of the time, while for CSB it is 47% and for 
RNS it is 99.1%. Thus, for gg, we usually expect the presence of four 6-jets in the final 
state (although some of these may fall below acceptance cuts or be merged with other 6-jets 
etc.). In the CSB case, the branching to t and 6 quarks is only mildly enhanced since all six 
squark flavors are extremely heavy. In addition, in the mSUGRA and CSB cases, gluinos 
only decay substantially to the lighter -ino states W\ and Z\_ 2 . For the RNS case, gluino 
decays to the light higgsino-like EWinos dominates but also decays to the heavier bino- and 
wino-likc states Z 3)4 and W 2 can be substantial. 

A diagram depicting gluino pair production followed by typical three-body decays is 
shown in Fig. [3j The presence of up to four 6-jets in the final state can be used as a 
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final state 

mSUGRA 

CSB 

RNS 

qq'Wi 

10.5 % 

34.0 % 

0.1 % 

tbWi 

13.4 % 

28.8 % 

45.6 % 

tbW 2 

-% 

-% 

2.2 % 

qqZ i 

3.1 % 

17.0 % 

-% 

bbZ\ 

0.5 % 

8.7 % 

-% 

ttZ\ 

60.3 % 

6.2 % 

17.2 % 

qqZ 2 

5.2 % 

2.4 % 

-% 

bbZ 2 

4.3 % 

0.3 % 

-% 

tiZ 2 

2.5 % 

3.0 % 

22.5 % 

tiz 3 

-% 

-% 

10.6 % 

ttz 3 

-% 

-% 

1.0 % 


TABLE II: Gluino branching fractions for the three benchmark models where q = u, d, c and s. 


powerful veto against dominant SM backgrounds such as tt production. Indeed, ATLAS 
searches[49] for gg production with > 3 6-jets in the final state offers the most powerful 
probe of gluino masses in the case where rrig -C rriq. 



C. Gluino cascade decay signatures 

We use Isajet 7.84[20] to generate a SUSY Les Houches Accord[50] (SLHA) hie for each 
benchmark scenario which is fed into Pvthia[5l] for generation of gluino pair production 
events followed by cascade decays. The gluino pair cross section is normalized to the NLO 
Prospino results of Fig. [2j We use the Snowmass SM background event set m for the 
background processes. The it background set is expected to be the dominant background[53]. 
where extra 6-jets can arise from initial/final state radiation and from jet mis-tags. While 
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the Snowmass background set was generated for y/s = 14 TeV LHC collisions, we have re¬ 
scaled the rates for y/s = 13 TeV collisions. Our signal and BG events are passed through 
the DelphesfM] toy detector simulation as set up for Snowmass analyses. 

We apply the following event selection cuts: 

• n(jets ) > 4 , 

• n(b-jets) > 3, 

• E T (j 1 ,j 2 - 4 ) > 100, 50 GeV, 

• for isolated leptons, then pr(£) > 20 GeV, 

• $ T > $ T [cut) = 50,100 — 500 GeV 

• A t > 1200 GeV, 

where A T = $ T + J2i ep tons E t + J2 je ts E t and for later use M eff = $ T + Eti E t(.U)- To 
gain some optimization of signal-to-background (S/B), we tried the above range of $ T cuts 
and evaluated S/B with and without the At cut. 

The cross sections after cuts for various multi-lepton + > 3 6-jets + $ T channels are 
shown in Fig. [ 4 } The optimal $ T cut for the Of* and 1 H channels was the hardest value: 
$ T > 500 GeV. For the Opposite Sign Same Flavor (OSSF) dilepton channel, the best cut 
was $ T > 400 GeV while for the Same Sign (S'S')-dilepton, 3H and 4H channels, the $ T > 50 
GeV was best. The A T > 1200 GeV cut helped just marginally. 

We see, from Fig. |dj that the signal cross sections after cuts in the jets+^ T (0H) channel 
are 1.9, 3.3 and 2.1 fb repectively for the mSUGRA, CSB and RNS cases while SM BG lies 
at 1.2 fb. In Fig. [5j we show the required value of LHC13 integrated luminosity which is 
needed to establish a 5cr signal, where in addition we also require at least 10 total signal 
events. From this plot, we see that just 8.3, 3.1 or 6.9 fb -1 of integrated luminosity L is 
needed to establish a first signal for the mSUGRA, CSB and RNS benchmark models with 
rrig = 1.4 TeV. The CSB benchmark model has a somewhat larger signal cross section and 
hence requires somewhat lower L in the OH channel as compared to the mSUGRA and RNS 
models since its decay modes include more hadronic and fewer leptonic cascades. 

In Fig. [4| we also see the cross section after cuts for the 1H channel. Even though one 
takes a leptonic branching fraction hit in this channel, the numerous sources for a single 
additional isolated lepton lead to cross sections after cuts which are comparable to those in 
the OH channel. For the 1H channel, RNS has the largest cross section 1.0 fb while mSUGRA 
and CSB are at the 0.8 fb level. This 1H + jets + $ T channel will confirm the signal which 
is already established in the OH channel with just a few additional (10-14) fb -1 of integrated 
luminosity. 

In Fig. [4] and [5j we also show the cross section after cuts and the required integrated 
luminosity for a 5cr signal for the OSSF, SS, 3H and 4H channels. These multi-lepton channels 
all exhibit a greater suppression due to multiple leptonic branching fractions as compared 
to the OH channel. For the 3H channel, background events come from isolated leptons in the 
6-quark decays. With the requirement of at least 3 6-jets, we do not observe any events in 
our ti background for the 4H channel. Also, in the multi-lepton channels, we see that the 
RNS model yields the largest cross sections due to the large gluino branching fractions into 
tops followed by t —* bW and W —> tv decay. From Fig. [HJ we see that typically ~ 100 fb -1 
is necessary to establish a signal in the dilepton and trilepton channels while rv_/ 10 3 - 10 4 
fb~ 4 would be required for a 5<r signal in the 4 1 channel. 
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FIG. 4: Cross section after cuts from gluino pair production for the three SUSY benchmark models 
and from tt background. For the 0 and 1 £ signals, we take $t > 500 GeV while for the OSSF 
dilepton channel we take > 400 GeV. For the SS, 3£ and 4£ signals, we require > 50 GeV. 



FIG. 5: Required integrated luminosity at LHC13 to establish a 5a SUSY discovery in various 
channels from gluino pair production for the three SUSY benchmark models compared to tt back¬ 
ground. For the 0 and If signals, we take > 500 GeV while for the OSSF dilepton channel we 
take $r > 400 GeV. For the SS, 3£ and A£ signals, we require > 50 GeV. 


IV. ESTABLISHING THE LSP ARCHETYPE 
A. Charged SUSY breaking 

One of the features of the CSB model is that W\ and the LSP are almost degenerate with 
a mass difference Am = — m^ i > m,±. W\ decays into charged pions at almost 100% 

rate. The reduced phase space also makes the chargino long-lived so that, once produced 
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FIG. 6: Distance to the decay point from the interaction vertex for the CSB benchmark model. 


at the interaction vertex, it travels a visible distance before it decays to soft pions plus 
the LSP. Since the chargino is so massive, its velocity is borderline relativistic leading to a 
highly-ionizing trail or track (HIT). The chargino lifetime is extracted from Isajet and 

the actual lifetime of each chargino is generated from the exponential decay law e ' w i. 
Then the track length is computed from d = /Lyf. 

In Figure |6j we display the histogram of the distance travelled from the interaction vertex 
to the decay point of each chargino. Here, we see that the typical length of each HIT is of 
order 2-20 cm. We also display the percentage of events containing 0-2 charginos. We see 
that 90% of the events passing our cuts contain either one or two charginos in each event. 
The presence of one-or-more HITs in candidate SUSY events would be the smoking gun 
signature of SUSY models with a wino-like LSP. 


B. Radiatively-driven naturalness (RNS) 

In the RNS benchmark model, it is emphasized[10) [19, 55) that the mass gap between 
the Z 2 and Z\ neutralinos is typically small: ~ 10 — 30 GeV which gives the inter-higgsino 
splitting. For our benchmark case, the value is Am = = 27.4 GeV. Notice 

this mass gap never gets much below about 10 GeV since naturalness also provides upper 
bounds to the gaugino masses via loop effects so that the higgsino-gaugino mass gap cannot 
become arbitrarily large. The modest Z 2 — Z\ mass gap has important consequences for 
phenomenology. It means that the Z 2 always decays via 3-body modes Z 2 —>■ Ziff which 
is dominated by Z* exchange. The decay mode Z 2 —* Z\£ + t~ occurs at 3-4% per lepton 
species, but the OSSF dilepton pair which emerges from this decay always has invariant 
mass kinematically bounded by — mg . This mass edge should be apparent in gluino 
pair cascade decay events which contain an OSSF dilepton pair. 

In Fig. [7J we show the invariant mass distribution of OSSF dilepton pairs in gluino pair 
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FIG. 7: Invariant mass of OSSF leptons. The dilepton mass edge and the Z peak are visible for 
the RNS model. We require n(b-jets ) > 3. 


cascade decay events where we require the above cuts but with $ T > max(100 GeV, 0.2M e //) 
and At > 1200 GeV and the presence of an isolated OSSF dilepton pair. The black his¬ 
togram shows the expected continuum background distribution arising mainly from ti pro¬ 
duction while the green histogram shows signal plus BG for the RNS benchmark model. 
The RNS signal is characterized by the distinct mass bump and edge below about 30 GeV. 
This feature provides the smoking gun signature for SUSY models with light higgsinosjlOj. 
One can also see a peak at m(£ + £~) ~ mz which arises from W 2 and Z 3]4 two-body decays 
to a real Z. The area under the m{£ + £~) < 30 GeV portion is ~ 0.025 fb so that of order 
400 fb -1 of integrated luminosity will be required before this feature begins to take shape 
in real data. 

For comparison, in Fig. [8] we show the same m(£ + £~) distribution for the case of the 
CSB benchmark. In the CSB case, first of all there are far fewer £ + £~ pairs present above 
background, and second there is no obvious structure to the signal distribution: we expect 
just a continuum. 

The second smoking gun signature for models with a higgsino LSP is the presence of same- 
sign diboson (SSdB) events which are from wino pair production [I9l 56] . In this case, the 
production reaction is typically pp —>- W 2 Z A followed by W 2 —> WiZ li2 and Z 4 —» W 4 W T . 
The Majorana nature of the Z 4 leads to equal amounts of same-sign and opposite sign 
dilepton events. Note that these SSdB events contain minimal jet activity- only that arising 
from initial state QCD radiation- as opposed to SS dilepton events from gluino and squark 
cascade decays which should be rich in the presence of additional high p T jets. 
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FIG. 8: Invariant mass of OSSF dileptons for the CSB model. We require n(b-jets) > 3. 
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FIG. 9: Invariant mass of OSSF dileptons for the mSUGRA model. We require n(b-jets) > 3. 


C. mSUGRA/CMSSM 

For the mSUGRA/CMSSM benchmark model with a 1.4 TeV gluino, then we expect the 
production of the usual multi-lepton+multi-jet + $ T cascade decay signatures as shown 
in Fig. [4j For the case of the mSUGRA benchmark, the mass gap between the wino-like 
Z 2 and the bino-like Z\ is 225.5 GeV so that the Z 2 -A Z\h (spoiler) decay mode is open. 
This two-body decay dominates the Z 2 branching fraction and so we expect no additional 
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structure in the dilepton invariant mass distribution. The m(£ + £~) distribution for the 
mSUGRA benchmark point is shown in Fig. [9] While no characteristic dilepton structure 
is apparent, it may be possible instead to pull out the presence of h —> bb decays in the 
mSUGRA cascade decay events where m(bb ) ~ m/57]. 


V. CONCLUSIONS 

During run 1 of the LHC at a/s = 7 — 8 TeV, the Standard Model was vigorously 
confirmed in both the electroweak and QCD sectors and the Higgs boson was discovered 
at rrih ~ 125 GeV. The presence of a bonahde fundamental scalar particle cries out for a 
mass stabilization mechanism of which the simplest and most elegant one is supersymmetry. 
Unfortunately, no SUSY particles have yet appeared leading to mass limits for the gluino 
particle of rrig >1.3 TeV. 

LHC Run 2 with a/s = 13 TeV has begun! New vistas in SUSY parameter space are open 
for exploration. While naturalness allows for gluinos as high as 4-5 TeV (with A ew < 30), 
it is yet true that naturalness (mildly via higher order contributions) prefers gluinos as light 
as possible. Motivated by these circumstances, we considered how SUSY discovery would 
unfold in three SUSY archetype models with a bino-, wino- and higgsino-like LSP each with 
a 1.4 TeV gluino, just beyond present bounds. 

We find that SUSY discovery could already arise at the 5cr level with just 3 — 8 fb -1 of 
integrated luminosity via the > 3 b — jet + $ T channel. Confirmation would soon follow 
in the >35 — jet + 1 — £ + channel. Further confirmation in the 2-3 lepton channels 
will require ~ 100 fb -1 . The CSB benchmark case would immediately be identified by the 
presence of one or more highly ionizing tracks in each signal event due to long-lived wino- 
like charginos which undergo delayed decays to a wino-likc LSP. No such HITs should be 
apparent in signal events from the mSUGRA or RNS archetype models. Instead, the RNS 
archtype would be signalled by a gradual buildup of structure in the m{£ + £~) OSSF dilepton 
mass distribution, where the m{£ + £~ ) < m^ 2 — mass edge along with a Z peak should 
be apparent with ~ 100 — 1000 fb -1 . In the RNS case, the gluino cascade decay events 
should ultimately be accompanied by the presence of same-sign diboson events arising from 
wino pair production. 

For the mSUGRA archtype with a bino-likc LSP, then we expect the usual assortment 
of gluino-pair-initiated cascade decay multilcpton+jets + $ T events but without HITs and 
without any apparent structure in the m(£ + £~) distribution. However, the presence of Higgs 
bosons lurking within the cascade decay events may be a distinguishing feature. On to data 
from LHC13! 
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